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Abstract. The theory of radiatively- induced warps in ac- 
cretion discs is applied to the discs of Be stars. It is found 
that these discs may develop warps in their inner regions, 
although once the warp amplitude is large enough then the 
interaction between the disc and fast radiatively-driven 
wind will determine its evolution. The warping is shown to 
be more important for later than earlier B stars. Although 
the interaction of the fast-wind with the disc will limit the 
amplitude of the warp, it cannot drive the warp radially 
outwards, and so the radial evolution of the warp depends 
on the dominant advective process within the disc. 

Typical timescales associated with growing modes are 
shown to be short, of the order of days-weeks, although 
these are not likely to be the timescales inferred from ob- 
servations of line-profile variations which are much longer, 
of the order of years. 



Key words: stars: emission-line, Be 
circumstellar matter 



stars: rotation 



1. Introduction 

Be stars are accepted to have a dense, slowly expanding 
disc in their equatorial planes and a fast radiatively-driven 
wind over their polar regions (see Slettebak 1988 for a re- 
view) . The emission line profiles generated in the disc have 
been monitored for many stars and are found to vary (e.g. 
Dachs 1987). Long-term variations in the ratio of the vi- 
olet to red components of a line, V/R, are attributed to 
m = 1 density waves in a Keplerian disc (Okazaki 1991, 
Papaloizou et al. 1992). Hanucshik (1996) also concludes 
that the disc is rotationally supported. Using the assump- 
tion that the disc is rotating at Keplerian velocities around 
the central star and is present throughout most of the 
star's life, Porter (1998) has recently showed that the ra- 
dial velocity at the inner edge of the disc must be small 
(v r <^ 0.01 km s^ 1 ). The conclusion of these studies is that 



the disc material orbits many times before it travels a sig- 
nificant radial distance - similar conditions to the material 
in an accretion disc. 

Pringle (1996), Malony, Bagelman & Pringle (1996) 
and Armitage & Pringle (1997) have recently examined 
the possibility that a warp may be induced in an accretion 
disc due to the incident radiation field from the central ob- 
ject. They find that if the accretion disc is optically thick 
at the wavelength at which the disc emission is maximal 
then the disc may become unstable. 

Here it is recognised that the physical situation exam- 
ined by Pringle and co-workers and that of the Be star 
disc are similar: the Be star disc is being illuminated from 
both sides by the central B star; the Be star disc mate- 
rial orbits many times before moving a signifiant radial 
distance, and observational as well as theoretical studies 
have shown that at some radius the Be star discs are opti- 
cally thick (e.g. Chokshi & Cohen 1987, Kastner & Maz- 
zali 1989, Marlborough et al. 1997). Recently, two Be stars 
have been examined and have shown observational charac- 
tertistics of a warped disc (Hummel, 1998). In this paper it 
is attempted to apply instability criteria for radiatively- 
driven warping to the discs thought to exist around Be 
stars. 

In §2 the disc structure is examined and the opti- 
cally thick regions are identified. The criterion for the 
radiatively-driven warp instability to arise in the disc is 
then examined in §3. The evolution of a warp interacting 
with fast radiatively-driven wind is considered in §4. This 
is discussed in §5, and conclusions given in §6. 

2. Disc structure 

The discs around Be stars have often been modelled with 
the "disc model" due to Waters (1986). This assumes that 
the disc density follows a power law in radius and, from 
mass conservation within the disc, that the radial velocity 
also has a power law form. Typically, the disc has a con- 
stant opening angle 8. Although this model is emperical, 
is is very successful in reproducing many of the observa- 



2 



John M. Porter: On the warping of Be star discs 



tions of Be star discs, (e.g. Waters 1986, Waters, Cote & 
Lamers 1987) 

However, theoretical models of disc formation and evo- 
lution have been less successful. Models have been devel- 
oped based on magnetism (Poe & Friend 1986), chang- 
ing wind line-driving parameters (Chen, Marlborough & 
Waters 1992), wind compression (Bjorkman & Cassinelli 
1993, Porter 1997), pulsation (Willson 1986) and viscous 
decretion (Lee, Saio & Osaki 1991). None of these mod- 
els may simultaneously account for many of the observed 
features of discs. Porter's (1998) small upper limit on the 
radial velocity of the disc at the star-disc boundary ap- 
pears to favour the viscous decretion disc model, but does 
not rule out other theories. 

Whatever the mechanism which produces the disc, 
its properties must be very similar to the disc model of 
Waters (1986); with the new limit on the radial veloc- 
ity (Porter 1998), the mass- loss rate in the disc is similar 
to that in the fast wind over the polar regions (also see 
Okazaki 1997). It is noted that as long as the viscosity is 
suitably defined then the disc model may be described us- 
ing identical mathematics as accretion discs (e.g. Pringlc 
1981) - indeed for the viscous decretion discs (Lee et al. 
1991) a viscosity is introduced from the outset. The main 
difference between the accretion discs and the decretion 
discs here (viscous or not) is essentially the boundary con- 
ditions - here the material drifts outwards; v r > and 
angular momentum is supplied from the star, possibly by 
pulsations (see Osaki 1986). 



2.1. Optical depth of the disc 



The model of Waters (1986) is used to represent disc 
structure. It is assumed to have a constant opening an- 
gle 9, and has a density profile of p = 10 -11 po(r / i?» )~™g 
cm~ 3 , where po is the density at the inner edge measured 
in KU n g cm~ 3 . The disc temperature is assumed to be 
Td = 0.5T o ff. It is also assumed to be isothermal, although 
this approximation may be easily relaxed (see Waters 1986 
for a discussion on non-isothermality) . 



at a radius R is calculated from the prescription in Waters 
(1986). Using his notation 



t v {R) = X x X td C R- 2n+1 



X X = \ 2 



hv 



(9ff + 9fb) 



X* d = 4.923 x 10 13 



7^3/2 2 



ipI 
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/ cos 2 "- 2 
J o 



V dy 



(1) 



where A is the wavelength in cm, gff and gjb are the Gaunt 
factors for free-free and free-bound emission respectively, 
z 2 is the mean squared atomic charge, 7 is the ratio of 
the number of electrons to ions, and p is the mean atomic 
mass. The opening angle of the disc is 9 and the rest of 
the symbols take their usual meanings. 

Using a selection of model stars spanning the main- 
sequence B star range various quantities for the model can 
be calculated (Table 1). The wavelength A m at which vF v 
is maximal is displayed in column 7 of Table 1 . To calcu- 
late the optical depth, it has been assumed that the base 
density in the discs for all the stars is identical: po = 1. 
There seems to be little correlation between the spectral 
type of the star, and disc tracers such as equivelent width 
of Ha (see e.g. van Kcrkwijk et al. 1995), and so the as- 
sumption of a common density normalization for all the 
stars is unlikely to be seriously in error. The calculation 
has also assumed that n — 2.5, and p = 1.6 -1 . 

Although the disc is optically thick in line transitions 
out to several tens of stellar radii (e.g. Chokshi & Cohen 
1987, Kastner & Mazzali 1989), eq.l does not yield r > 
1 for A m . Therefore the wavelength at which the optical 
depth is greater than unity A T =i at the inner edge of the 
disc is calculated, as is the fraction of emission at that 
wavelength to the maximal emission vF v . These quantities 
are presented in columns 8 and 9 of Table 1. The fraction 
of the disc emission at A r= i is small for the earliest B 
star, rising to 60% for the latest spectral type considered. 
The effect of this partial optical thickness of the disc is 
likely to be on the magnitude of the warp, should the 
disc become unstable (see below). The discs around the 
earliest (latest) B stars will generate the smallest (largest) 
amplitude warp. Clearly, if the disc is to become unstable 
to radiation-induced warping, it must be unstable to these 
perturbations very close to the star. 



From this the optical depth of the disc may be calcu- 
lated at a frequency v at which vF v at the disc temper- 
ature Td is maximum (where F v is the Planck function). 
Using cylindrical geometry, the optical depth of the disc r„ 



3. Disc warping 

The ciriterion that radiation-driven warping has been de- 
rived by Pringle (1996) and Malony, Bagelman & Pringlc 
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Table 1. Stellar parameters for the stars considered from Schmidt-Kaler (1982). The mass loss rates (in Mq yr~ ) come from 
the relation by Garmany et al. (1981). A m is the wavelength at which vF v at the temperature of the disc (T d — 0.5T c g) is 
maximal. The wavelength at which the disc has an optical of unity at the inner edge is A T =i. The fraction of disc emission at 
A T= i to that at A m is listed in the ninth column. The tenth column gives the growth timescale ir for the warping instability at 
r = R* from eq.7. The final column gives the bending (from eq. 12) at which the wind-disc interaction becomes significant. 



Spectral 


M» 




logL* 




logM„ 


A m 


A r = l 


emission 


tr 


ip-6 


type 


(M e ) 


(Re) 




(10 4 K) 




(H 


(Mm) 


fraction 


(days) 


(degrees 


BO 


17.5 


7.4 


4.81 


33.9 


-7.5 


0.22 


1.29 
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2.45 
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-11.6 
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1.22 
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20 
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3.4 
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2.28 


11.2 


-11.8 


0.67 


1.18 


0.60 


8.0 


18 



(1996). The typical timescale for the radiation torque tr 
and the viscous timescale for damping warps t U2 are: 



t U2 = —, t r = 12ttG 1/2 c 

V2 



(2) 



where v-i is the (R, z) component of the disc viscosity (see 
e.g. Pringle 1992) M* and L* are the mass and luminos- 
ity of the star, G and c are the gravitational constant 
and speed of light respectively, and S is the disc sur- 
face density. For instability, Pringle (1996) estimated that 
> 27rfr, a result which has been found consistent with 
numerical simulations (Pringle 1997, Armitage & Pringle 
1997). 

To make further progress on this problem for decre- 
tion type discs, however, the angular momentum transport 
equation must be examined again. For a steady disc, a con- 
stant C is produced when the angular momentum equation 
is integrated (Pringle's 1981 eq.3.8). For decretion discs, 
this constant is the rate of angular momentum input at the 
inner boundary of the disc. If the central star rotates with 
an angular velocity of a fraction /* of the Keplerian value, 
then the constant C should become C = (1 — f*)MR 2 Q. 
This leads to a modification of Pringle's (1981) eq. 3.9: 



M 

v{L = — 
3tt 



(1 



1/2 



(3) 



where vi is the (R, <j)) component of viscosity. 

With this and the instability criterion eq.2 then be- 
comes 



R \ 1/2 
it) 



LR\' 2 



knG x l 2 cr\M^I 2 M 



<0. (4) 



This is essentially different from the expression of Ar- 
mitage & Pringle (1997) (their eq.3) as the full expression 
of our eq. 3 has been used instead of the approximation in 



the case of large R/R*. As the disc is (partially) optically 
thick only in the inner parts then the full expression is 
required. 

Setting eq.4 equal to zero then defines the critical 
radius R c separating regions of instability and stability. 
With the replacement r = R c /R* , this becomes 



„l/2 



3.7 x 1CT 2 



(l-h)-gr = Q 



r)M. 



10 



(u \ (m* \ _l/2 ( r* y/ 2 



(5) 



where T] = 1/2/vi, and the mass loss rate is measured 

in lO- lo M yr- 1 . 

This equation yields Armitage & Pringle's result (their 
eq.5) if the (1 — /*) term is absent. Eq.5 is a quadratic 
equation for r 1 / 2 , and so it may produce zero, one or two 
solutions depending on the value of g. The solution is 



1 ± y/1 - 4ff(l - gj 
25 



(6) 



If only one solution is produced, then it is the inner bound- 
ary beyond which warping instabilities may occur; if two 
solutions are produced (both with r > 1), then the re- 
gion between them is stable to the warping instability. If 
4g(l — /*) > 1 then no real solutions arc obtained. This 
occurs when eq.4 is obeyed for all values of r, i.e. the whole 
disc may be unstable to warping - in this case the warping 
occurs everywhere the disc is optically thick. 

To calculate the function g the mass loss rate in the 
disc and 77 needs to be supplied. This is problematical for 
both quantities. As the viscosity may not be isotropic (see 
Pringle 1992), r\ may not necessarily be unity. Also, al- 
though the disc density may be derived via observations 
of Be star disc continuum emission in the IR (see e.g. 
Waters 1986) reasonably accurately, the radial velocity in 
the disc, and hence the mass-loss rate is unknown. De- 
cretion disc models produce low radial velocities (see e.g. 
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Okazaki 1997), and from angular momentum considera- 
tions the radial velocity at the stellar surface must be 
v r <; 0.01 km s -1 (Porter 1998, but see Porter's discus- 
sion). Low velocities such as these imply mass- loss rates 
in Be star discs similar to those for the radiation-driven 
wind over the poles (Okazaki 1997). However, if the disc 
is a viscous decretion disc, then there is no a priori rea- 
son that its mass loss rate and the fast wind's mass loss 
rate should be the same as two different mechanisms are 
acting. 

Clearly, there is some uncertainty here. To keep the 
calculation general, an "ignorance" parameter A is intro- 
duced such that nM_ w is replaced by AriM w , where M w 
is the mass-loss rate of the fast wind over the pole (taken 
to be that from Garmany et al.'s 1981 emperical relation; 
column 6 Table 1). With this introduction, the function 
4g(l — /„) is calculated for the sample of Be stars and is 
shown in fig.l. 
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/1000 







BO Bl B2 B3 B4 B5 B6 
spectral type 



B7 B8 



Fig. 1. The function 4g(l — /,) defined in the text. The four 
lines are labelled with the different values of A they correspond 
to. The dotted line is the limit above which no real solutions 
are obtained for eq.6. 



Fig.l shows that unless A is large, the disc is unstable 
for all radii where it is optically thick. If the mass loss rate 
in the disc is large, or the viscosity is very anisotropic, then 
eq.6 does have two real solutions, and a region of the disc 
will be stable. It is found that the inner solution, where 
it exists, is located within the star (r < 1), and therefore 
only the second solution is meaningful. This outer solution 



then provides a minimum radius at which the disc becomes 
unstable. As an example calculation, the outer solutions 
are calculated, where they exist and r > 1 for values of 
A = 100, 1000 and are displayed in fig.2. 




B0 Bl B2 B3 B4 B5 B6 

spectral type 



Fig. 2. The outer solution to eq.5, for the labelled values of 
A. There are no solutions which have r > 1 for A < 12 for the 
stellar parameters in Table 1. 



3.1. Growth and precession rates 

How fast do the unstable modes grow when they are ex- 
cited? The growth timescale from eq.2 tr in years is 



t r = 1.4 



R* 

R Q 



2 fM* 



Ma 



po tan# 



R 
R* 



(7) 



The precession timescale of the warp is T w Antr which, 
when R/R* « 1, po w 1 and tan# « 0.1 becomes 



T w 1.7 



Rq 



yr- 



(8) 



The growth timescales are displayed in column 10 of Ta- 
ble 1. It should be pointed out that these timescales 
are short, typically O(days). Consequently, prcccsion 
timescales are O(months). The observational implications 
of these timescales is returned to later in the discussion. 
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4. Disc fast wind interaction 

Once the disc starts to warp, it will interact with the fast, 
radiatively driven wind from the star (which has radial 
velocities of ~ 10 3 km s _1 , see e.g. Snow 1981). A bending 
mode will present a non-negligible area to the fast wind 
and the fast wind will interact via its ram-pressure. In 
order to attempt to describe the evolution of the warp in 
the optically-thin regions of the disc where the radiatively- 
drivcn warping mechanism does not operate, some typical 
timescales are now calculated. 

4-1. Latitudinal timescale and magnitude of the warp 

The fast wind will interact with the warped disc and at- 
tempt to deflect the warp toward the plane of neutral sta- 
bility - the equatorial plane. This interaction is complex, 
however, as it depends critically on the amplitude, (or lo- 
cal inclination) of the warp, and is therefore a function 
of the warp itself. The timescale for this deflection is de- 
noted as the wind-disc timescale t w d and is now estimated. 
The local disc opening angle is ip, (for bends ip is a mono- 
tonically increasing function of R), and the opening angle 
of a stationary disc is 9 (see Fig.l, Porter (1997). The 
velocity of the wind perpendicular to the disc surface is 
v± = v w sm(tp — 6), where v w is the velocity of the fast 
wind (which is assumed not to have a meridional com- 
ponent - although see Bjorkamn & Cassinelli 1993, and 
Owocki, Cramner & Gayley 1996) The ram pressure di- 
rected toward the equatorial plane is then p w v \ where p w 
is the density of the wind. This pressure acts on an area 
of $i? SR/ costp, where $ is the angle over which the warp 
exists (e.g. for a simple inclined disc, <f> = w). 

The acceleration g w d toward the equatorial plane pro- 
vided by this interaction is equal to the force due to the 
ram-pressure (resolved into the z direction), divided by 
the mass dm on which it acts dm — <&R SRpdH, where pd 
is the disc density and H = R tan# is the scale height of 
the disc. Hence, the acceleration is: 



p w vl \ sin 2 (V> - 6) 



Pd 



Rt&n9 



(9) 



The timescale for "flattening out" warps by this wind- 
disc interaction depends on how far the disc has moved 
from the equatorial plane 5H = R(t&nip — tan#). Assuming 
that both and ip are small, then the timescale is t? wd ~ 
R{i> - 6)/9wd, or 



^wd R 



Pd 

PwVl 



9 



ip-6 



(10) 



Setting p u 



M w /4ttR 2 v w , and again using pd 



10 11 po(R/R^ ) n , this timescale (in seconds) becomes 



t wd = 6.8 x 10 4 



R_ 



2-# 



\R Q 



P6 



vIm: 



(ii) 



where v$ is the wind velocity in 10 8 cm s 1 (10 3 km s 1 ) 
and the wind mass-loss rate M w is measured in 
lO^MQyr- 1 . 

This expression may now be used to estimate the mag- 
nitude of bending modes which may grow via the effect of 
radiation before they are seriously peturbed by the wind- 
disc interaction. This is done by calculating the bending 
(ip — 9) at which the wind-disc timescale t w d is equal to 
the radiatively-driven warping timescale tr . This yields 



6{ip -9) = 1.6 x 10" 



Mr., 



( R* \ I L 



1 



p v s M v 



R Q 
R 

R* 



n-1 



(12) 



Assuming «g«l, and values of M w from Garmany et al.'s 
(1981) relation, the values of ip — 9 arc listed in column 9 
of Table 1 at the inner edges of the disc (-R/-R* ~ 1), for 
9 = O.lrad - it is clear that large warps can develop before 
they are seriously affected by the wind-disc interaction. In 
fact for the early B stars the approximation that ip is small 
is severely strained. 

4-2. Is the warp driven to larger radii 

Clearly, the fast wind will attempt to force the warp to 
larger radii. The effectiveness of this depends on the ratio 
of the radial momentum flux in the wind to the gravita- 
tional attraction of the disc. If this ratio is greater (less) 
than unity then the warp will (will not) be driven out- 
wards. 

The wind's momentum flux is simply the product of 
mass-loss rate and the wind velocity, weighted by the ratio 
of the solid angle of interaction to 47rsteradians. The area 
of interaction is &RSRsm((p — 9)/(cos9 cosi/0, and so the 
wind's momentum flux T is 



= M / SRsH^-0) 
\2R cos9 cos?/; 



(13) 



The gravitational interaction of the star and the mass over 
which the wind's momentum flux acts is 



_ GM» $GAf„ p d HSR 
T g = —d m = . 



(14) 



Scaling the wind's mass-loss rate and velocity to 
10~ 10 M Q yr _1 , and 10 8 cm s _1 respectively, and using the 
power law expression for the disc, the ratio of the two 
terms above is 



^ = 5.4 x 10- 4 M w v 8 

■J a 



Mr, 



R* \ Y_R 

sin(V> — 9) 
sin0 cos^ 



(15) 



In the last section it was found that ip will be tens of 
degrees when the wind-disc interaction takes placeand so 
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the trigonometric term will numerically be ~10s. There- 
fore, it is clear that the fast wind cannot drive the warps 
to large radii. 

5. Discussion 

From the preceding sections it seems that Be star discs 
do posess the correct attributes to enable them to become 
unstable to radiatively-driven warping. If this is so, then 
they should be observable. In fact it may well be that 
this is indeed the case - Hummel (1998) has indeed inter- 
preted line-profile variations in 7Cas (BOIVe) and 59 Cyg 
(B1.5Vnne) to be due to a warped inner disc. This adds a 
large degree of credibility to the above approach. Hummel 
did consider the radiatively-induced warping mechanism 
as an explanation for the observations, but rejected the 
mechanism as he found a prohibitively large mimimum 
radius for the instability to operate (10 7 i?* ). The incon- 
sistency between his result and that presented here lies in 
the application of accretion disc theory to decretion discs 
of Be stars. (This involved approximating the accretion 

efficiency e = L/Mc 2 by the ratio of the Schwarzschild ra- 
dius to the stellar radius, leading to an error of ~ 5 x 10 3 . 
This value is squared to yield the inner radius for insta- 
bility, which gives an overestimate of ~ 2 x 10 7 , bringing 
Hummel's value into line with that here.) 

A feature of the development of the warps which is 
specific to Be stars is the action of the central star's ra- 
diatively driven wind. This introduces an extra (warp de- 
pendant) effective viscosity leading to a different evolution 
of bending modes. 

5.1. Development of a warp 

In the region of the disc where it is optically thick, ra- 
diatively induced warps grow on the timescale tr- These 
warps are unmolested until they become large amplitude 
bends (defined from eq.12). At this point the wind-disc in- 
teraction will start to domiate over the induced warping, 
and the bending mode will oscillate about the equatorial 
plane under the action of the fast wind's ram pressure. 
If it is still within the optically thick part of the disc, 
then the combined effects of the induced warping, and the 
wind-disc interaction will limit the oscillation's amplitude. 
However, when the mode propagates out of the optically 
thick part of the disc, then the radiation-induced warping 
becomes inoperable, and the evolution of the mode is de- 
termined entirely by the wind-disc interaction. It has been 
found that the warp will not be driven to large radii un- 
der this interaction, and will only propagate outwards via 
the action of viscosity (for decretion discs) or advection. 
However, it is difficult to determine whether the warp will 
be amplified or whether it will decay due to the effects of 
the wind - further work on this point is necessary (Porter, 
in preparation). 



5.2. Observational signatures 

The timescales for the growth of radiatively-driven warps 
are shown in Table 1. Although modes may be excited 
over this short timescale, it is not necessarily the timescale 
over which observed disc tracers indicate the warp. Once a 
mode propagates outward beyond the optically thick part 
of the disc, then the timescale for its evolution and preces- 
sion will be determined by the wind-disc timescale (eq.ll). 
Consider a bending mode grows to its maximum ampli- 
tude where the wind-disc interaction becomes effective. 
Once it has propagated out in the disc and its amplitude 
has decreases, then the timescale the mode evolves at falls 
(see eq.ll). 

The observational properties of such a mode are not 
defined simply by the optically thick part of the disc de- 
fined in §2. Typically, the hydrogen lines form within sev- 
eral tens of stellar radii (e.g. Hanuschik et al. 1988) and 
it is these which are commonly used to infer properties 
about the disc. Consequently, the warps inferred from the 
observations are likely to have longer timescales than the 
growth timescales, and precession rates in eq.7 and 8. 

For example, consider a warp which had propagated 
to 10i?* , and has ip — 9 = 0.1. The precession timescale 
for the wind-disc interaction T — 4Trt wc [, ranges from 170 
days for a BO star to 4200 days for a B8 star, consider- 
ably longer than the timescale over which the mode grows. 
Therefore it is likely that observationally derived preces- 
sion timescales will range from several hundreds, to several 
thousand days. This is broadly consistent with Hummel 
(1998), who ascribes timescales of <~ lOOOdays for preces- 
sion of the discs around 7 Cas and 59 Cyg. 

6. Conclusion 

The theory of radiatively-induced warping of discs has 
been applied to Be stars and found to be important, espe- 
cially for late B stars. The timescales associated with the 
growth of warps is short O(days), although these may not 
be the relevant timescales over which the modes are ex- 
amined observationally. It is noted that warping has been 
observed in Be star discs already (Hummel 1998), and 
that the observationally derived precession timescales fall 
in the range in which the simple argument predicts. 

A scenario has been presented in which radiatively- 
driven warps are created in the central parts of the disc 
and propagate outwards under the influence of the inter- 
action of the disc and fast stellar wind. 
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